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Abstract 

If the lepton numbers are associated with global symmetries spontaneously 
broken at a scale below 1 TeV, neutrino oscillations in supernovae produce 
classic Majoron fields that perturb the neutrino propagation itself and may 
change the oscillation patterns in the periods of largest v fluxes. The impact 
of the Majoron fields on the same transitions as v e — > ux that presumably 
occur in the Sun is studied in the case of the non-adiabatic MSW solu- 
tion. It is shown how the back reaction of the Majoron fields may improve 
the adiabaticity of these oscillations in a supernova environment which has 
implications on the outgoing u e spectrum. 

PACS numbers: 14.80.Mz, 14.60.Pq 
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I. INTRODUCTION 



Several experiments from solar to atmospheric neutrinos and laboratory oscillation 
experiments |1] indicate that neutrinos oscillate and leptonic flavors are not conserved. 
This is in contrast with the minimal standard model (SM). It is quite possible that the 
energy scale of breaking of the three lepton numbers is comparable or even smaller than 
the Fermi scale. Furthermore, if they are spontaneously broken by the expectation values 
of some scalar fields then, a few bosons with zero mass should exist - the Nambu-Goldstone 
bosons (NG) - one per broken global symmetry. 

It was recently pointed out that these NG bosons (called Majorons or familons 
when associated with lepton numbers) couple to the time rate of creation of the respective 
lepton numbers carried by the matter particles and therefore coherent NG fields are 
produced whenever lepton number violating processes occur simultaneously. That is the 
case if neutrinos change flavor on their way out from stars as seems to happen in the 
solar system. Once the NG fields are generated (the triggering process may be a normal 
Mikheyev-Smirnov-Wolfenstein (MSW) resonant conversion Q), they change in turn the 
relative potentials of the different neutrino species and so get a life of their own. 

The numbers show that if the scale of symmetry breaking is below 1 TeV then, the 
Majoron fields are important enough to play a role in supernova neutrino oscillations. 
They are however too small in the case of the Sun unless the scale of symmetry breaking 
lies below 1 KeV. As a result, supernova neutrinos may exhibit oscillation patterns in 
contradiction with the observations of solar, atmospheric and terrestrial neutrinos. We 
ought to be prepared, in the event of a close by supernova explosion, for the possible kind 
of effects caused by NG fields. 

In the previous paper f2|, the example studied was that Majoron fields generated by 
the conversion v e — > u T , assumed to take place in a certain resonance shell, yield neutrino 
potentials which become competitive with the standard electroweak potentials at larger 
radii and therefore affect the other flavor transitions characterized by smaller Am 2 . The 
effects can be so dramatic as the resonant oscillation v e <-> in a context of m Ve < m Ufi 
hierarchy, where the resonance is otherwise possible for v e <-> z/ M but not for the anti- 
neutrinos, if they only interact via standard W and Z bosons. In the present paper, I 
want to discuss what happens if the Majoron potentials are already significant in the very 
region where the oscillations they are generated from occur. Then, a back reaction effect 
takes place yielding an interesting flavor dynamics. 

Suppose that in the Sun the electron-neutrino oscillates into the muon-neutrino with 
the parameters of the non-adiabatic, small mixing angle solution [f|||. Then, the v e <-> 
transitions are also non-adiabatic in a supernova and only a fraction of each neutrino 
species is converted into the other. Furthermore, since the level crossing probability is an 
increasing function of the energy, the hotter z/ M s have larger survival probabilities than 
the cooler u e s. It will be shown that the back reaction of the NG fields improves the 
adiabaticity of the neutrino transitions, thus yielding a hotter energy spectrum for the 
outgoing u e s. That kind of effect can in principle be traced in those detectors such as 
Super-Kamiokande and SNO, capable of detecting supernova electron-neutrinos ||. 
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II. MAJORON AND NEUTRINO FLAVOR DYNAMICS 



In the following it is assumed that the partial lepton number L e is conserved at the 
Lagrangian level but the global symmetry associated with it, U(l) e , is spontaneously 
broken by the expectation values of one or more scalar iso-singlets <7j. Then, a NG boson 
£ e exists with zero mass. The neutrino mass matrix violates in principle the three lepton 
numbers but for simplicity I will ignore the other possibly existing Majorons. It may 
be interpreted as meaning that the respective scales of symmetry breaking are slightly 
higher. 

It is well established that the Nambu-Goldstone bosons only interact through deriva- 
tive couplings |7]|| (related to the soft pion low energy theorems). This is clarified |].[^] 
by changing variables from the original fields with definite L e charges namely, fermions 
X a and scalars a,, as follows: 

X a = exp(-i£ e L a e )r , (1) 

a t = exp(-z£ e L l e ) ((a*) + . (2) 

The so defined physical weak eigenstates, the fermions ip a = e, u e , ... and massive bosons 
pi, are invariant under the group U(l) e . In terms of these fields the symmetry is simply 
realized as translations of the Majoron field, £ e — > £ e + a, and consequently the Lagrangian 
can only depend on the derivatives d^ e . 

The non-standard interactions relevant to this work are contained in the expression 

£ = - v f C m ab v\ + h.c. + ]pl d^ e d^ e + C iQt (d,Q , (3) 

where VL 2 e = 2J2i \L\ (cr^) | 2 and z/£ denote the standard v e Ll v T L or any extra neutrino 
singlets. The £ e equation of motion, 

d,d^ e = -d,Jt?/n 2 e , (4) 

identifies with the conservation law of the Noether current associated with the symmetry 
£e ~ *• + a. All the Majoron interactions are cast in the current J£ obtained from 
Ant(d M £ e )- Its leading terms do not depend on the particular model and are derived from 
the x° \ °i kinetic Lagrangian using Eqs. (HD, (Q). The result is the following: 

Ant = (d^ e ) (e7 f, e + v ( r 1 ii v e + • • •) , (5a) 

JM= eYe + ^ e + - • (5b) 

The dots represent scalar boson terms and model dependent radiative corrections that 
are not relevant for this work and will ignored in the following. 

The expressions above are typical of one kind of models, the Abelian-singlet Majorons 
|10|-|12],3 . Singlet means that the scalar fields that spontaneously break the lepton number 
symmetries are all singlets under the SM gauge group SU(2) xU(l), thus complying with 
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the LEP results on the Z° invisible width, unlike the triplet-Majoron model [13]. By 
Abelian I mean U(l) symmetry groups, not horizontal in flavor space, associated to L e , 
L^, L T or any linear combinations of them. As the respective currents do not change 



flavor these models are not bound by the laboratory limits on the familon models JT4 
In addition, of all interactions involving SM particles the neutrino masses are the most 
important source of lepton numbers violation. Thus, in single collision or decay reactions, 
the effective strength of the Majoron couplings, resulting from d^J^/ (o~) , is proportional 
to the neutrino masses and suppressed by the symmetry breaking scale: g ~ m„/ (a). For 
that reason, one does not expect observable zero neutrino double beta decays accompanied 
by Majoron emission as shown in ||15||. With a sensitivity to neutrino-Majoron couplings 



of the order of 10 5 0] they cannot even probe relatively low symetry breaking scales. 



This kind of models also escape present astrophysical bounds on the couplings of Nambu- 
Goldstone bosons. Neutrino-Majoron couplings with strengths g ~ m v j (a) are too far 
from 10~ 4 to change supernova collapse dynamics [JT7H , and even below the threshold of 
~ 10~ 8 ' 5 for supernova cooling through singlet Majoron emission [D|. Finally, the pseudo- 



scalar couplings to electrons, that could be responsible for energy loss in stars |T!|, only 
arise through radiative corrections and are so further suppressed. 

In the £ e equation of motion, the source term d^J£ is nothing but the time rate of 
creation of L e -number carried by matter particles per unity of volume. If the neutrinos 
v e and oscillate into each other outside the supernova neutrinospheres, but not their 
anti-particles, the net variation of L e is given by the difference between the numbers of 
converted neutrinos N(v^ — > u e ) and N(u e — > i/^). In a stationary regime the fluxes are 
constant in time and Eq. (|J) reduces to a Poisson equation with a Coulombian solution for 
£ e [[|. The gradient A e = — V£ e obeys a Gauss law. In a spherical symmetric configuration 
it only has a radial component, 

that is determined by the integral of the source term over the volume of radius r, in this 
case, L e -number created per unity of time in that volume, L e (r) = J d 3 x d^Jg. It can 
also be expressed as [Q] 

Ae ( r ) = ~~& -> Ve) - j("e -> ^)] , (7) 

where j(u e — > v^) denotes the flux of e- neutrinos converted to \x flavor and jiy^ —>■ u e ) the 
reciprocal. Both these quantities are functions of the radius r. 

Taking in consideration the interactions specified by Eqs. (|J), (|5aD , the equation of 
motion of the neutrino wave function is, including flavor space, 



7°V - i)L = mil) R 



where m is the v mass matrix (real for simplicity), L e is the flavor- valued quantum number 
(1 for u e and for v T ) and Vq designates the flavor conserving SM potential in a medium 
at rest. One derives in the same way as in the case of a scalar potential Vq the equations 
governing flavor oscillations |2IJ (see also 0), 
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+ V + A e L e 



(9) 



which give (after absorbing a flavor universal term in the wave function) 

/ 2E(V W + A e ) \Am 2 sin 26 \ I v P 

\ |Am 2 sin26> Am 2 cos 26 



.d_ 

dr 



V,- 



Vn 



1 

2E 



(10) 



V w = V2G F n e is the charged- current potential of v e in a medium with electron number 
density n e [^TJ and 6 is the mixing angle. It is worth to notice that the results do 
not change if one considers alternatively a NG boson associated with breaking of L M or 
L e — Lp. The reason is, v e <-> oscillations only violate L e — (the non-conservation of 
L e + L M is suppressed by m 2 /E 2 ) and only care about the difference between the u e and 
Up potentials. 

In order to calculate A e (r) one needs to know the fluxes of and v e as functions 
of the radius with and without oscillations. I make the simplification of neglecting the 
oscillation length by saying that a neutrino with energy Er is converted to the other flavor 
at the position where the resonance condition, 



E, 



Am 2 cos 26 
2{V W + A e ) ' 



is fulfilled. In addition, having in mind that in a non-adiabatic regime only a fraction 
1 — P c changes flavor (small mixing angle), the level crossing probability P c is calculated 
using the Landau-Zener approximation pL|2(J, 



P C {E) = exp 



7T 

4 



Am 2 sin 2 26 



dE R /dr cos 26 



'12) 



IR 



I believe that these approximations change the numbers but not the lesson taken by 
comparing the results with and without a Majoron field. 

Let the number of emitted particles per unity of time and energy be specified by the 
distribution functions fv e (E) = dN u JdE and f v (E) = dN v /dE. They are normalized 
by the relation between the number luminosity A^^, the energy luminosity L u and the 
average energy E v of each v flavor namely, N v = L u /E u . L u and N u will be given below in 



unities of ergs/s and ergs/s/MeV respectively. Equation (|TTD gives the resonance energy 
as a function of the radius: particles with lower energies reach the resonance position 
at higher density regions and the hottest neutrinos oscillate at the largest radii. The 
statement that the number of converted neutrinos is the fraction 1 — P c of the number 
of particles with resonance energy translates into an equation for the time rate dL e = 



dN(^ 



dN{v e 



dL f 



of L e -number creation in a shell with depth dr: 



Pr 



L (Er) - fv e {Er] 



dEj 
dr 



-dr . 



(13) 



This establishes a differential equation for L e (r). Notice that the derivative of Er is not 
independent from the derivative of L e because the Majoron potential A e (r) that enters in 
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FIG. 1. Level crossing probability P c as a function of the v energy. Dashed curve holds 
for the SM with constants M = 4 x 10 31 g and Y e = 1/2. In the dotted and bold curves the 
Nambu-Goldstone field £ e exists with G e = and G e = 4 Gf respectively, and the luminosities 
are 10 52 ergs/s for v e and 7 x 10 51 ergs/s for u^. 



E R depends also on L e (r), as Eq. @ shows. In addition, the probability P c depends on 
the Er and L e derivatives as well, and that makes a non-linear differential equation for 
L e {r) specified by Eqs. ©, (0) - (0). 

It remains to tell the density profile of the medium. In the regions of a supernova star 
with densities typical of the Sun the mass density goes as 1/r 3 , the constant M = pr 3 
lying between 10 31 g and 15 x 10 31 g depending on the star [T2^]. Then, in terms of the 
electron abundance Y e w 1/2, the electroweak potential (V2Gpn e ) is 



V w = 0.76 Y f 



(14) 



where rio = r/10 cm. This is to be compared with 



A P = 1.48 



G F 10 51 ergs/s/MeV 



.-2 
10 



x lfr 12 eV 



(15) 



where Gf = 11.66 TeV -2 is the Fermi constant and G e = l/f^- Clearly, if the neutrino 
luminosities are sufficiently high say, 10 52 ergs/s for 10 MeV neutrinos, and the scale of 
lepton symmetry breaking is around or below the Fermi scale, the Majoron potential 
A e becomes competitive with Vw at radii where the resonance occurs for Am 2 values 
interesting for solar neutrinos (10~ 5 — 10~ 4 eV 2 ). More generally, at large enough distances 
the Majoron potentials decay as the inverse square radius and necessarily dominate over 
the local interactions. 
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FIG. 2. Total potential V e = V Ue — V v as a function of the radius. The dashed curve stands 
for the SM potential, the dotted and bold curves for the Majoron case with the same parameters 
as for the homologous curves of Fig. 1. 



Let us examine the v e <-» oscillations with the mixing parameters of the non- 
adiabatic solar neutrino solution (for a recent update see P^j), choosing in particular the 
values Am 2 = 7 x 10 _6 eV 2 , Am 2 sin 2 29 = 4 x 10~ 8 eV 2 . In a supernova the resonance 
is non-adiabatic as well and, as Eq. (0) indicates, the survival probability P c increases 
with the v energy. This was studied in detail in the framework of the SM | |24|| . In Fig. 1, 
P c in the Landau-Zener approximation is plotted against the energy. The dashed curve 
holds for the SM potential with a constant M = 4 x 10 31 g. It is manifest the aggravation 
of the non-adiabaticity with the energy. 

To study the Majoron case one has to specify the energy spectra and luminosities. 
I used Fermi-Dirac distributions with the following values of temperature and chemical 
potential Q: for u e , T = 2.4 MeV and fi = 3.2 T; for T = 5.1 MeV, \i = 4.1 T. This 
gives average energies of 10 and 23 MeV respectively. The luminosity intensities are in 
turn 10 52 ergs/ s for v e and 7 x 10 51 ergs/ s for which amount to particle emission rates of 
10 51 and 3 x 10 50 ergs/s/MeV respectively. Because the e-neutrinos are more numerous, 
the v e <-> i/p oscillations produce a net destruction of L e -number and a positive Majoron 
potential A e . The dynamics is the following: the less energetic u e oscillate to at the 
smallest radii; this conversion produces a positive A e (r) which attenuates the fall of the 
total potential V e = Vw + with the radius; as a consequence, the adiabaticity improves 
at larger r and the most energetic neutrinos change flavor with higher probabilities. 

Figure 2 shows the total potential 14 as a function of the radius. The dashed curve 
stands for Vw, the SM potential, with M = 4 x 10 31 g and Y e = 1/2. In the dotted and 
bold curves the Nambu-Goldstone field associated to L e symmetry breaking operates with 
a constant G e = Gp and G e = 4Gp ; respectively. The potential falls more slowly as the 
field £ e grows. The level crossing probability is plotted in Fig. 1 for both cases (dotted and 
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FIG. 3. The dotted curves are the assumed luminosity distributions for u e and as emitted 
from the neutrinospheres. The dashed curve represents the u e luminosity after electroweak 
neutrino oscillations whereas in the bold curve a £ e field exists with G e = 4Gf- 

bold curves) and the effect is clear: the stronger the Majoron field, the more efficient is the 
flavor conversion. It is worth to mention that if the /z-neutrinos were more numerous than 
the e-neutrinos the effect would be the opposite because the Majoron potential would be 
negative (L e > 0). That is actually reflected in the rapid rise of P c at the energy band 
around 20 MeV. 

Figure 3 shows the implications for the outgoing v e energy spectrum. The dotted 
curves are the assumed luminosity distributions for the emitted u e s and z/ M s. The dashed 
curve represents the luminosity of the e-neutrinos that come out of the star after standard 
MSW oscillations and the bold curve is the same but with a Majoron field (G e = 4G F )- 
The improvement of adiabaticity makes more z/ M s to convert into v e and less u e to survive, 
and because the /z-neutrinos are more energetic, the outgoing v e spectrum is harder than 
if there was no Majoron field. The average energy of the outgoing e-neutrinos is 13 MeV 
if G e = but rises to 17 MeV if G e = G F and 21 MeV if G e = 4 G F . 

III. CONCLUSIONS AND DISCUSSION 

To summarize, if the explanation of the solar neutrino deficit is the MSW non-adiabatic 
oscillation v e — > (or v e — > u T ) then, the standard model of electroweak interactions 
predicts that in a supernova the v e <-> transitions are also non-adiabatic. It means 
that, to a large extent, the e-neutrinos preserve their lower energy spectrum, unless v e also 
mixes to another flavor with a too high or too low Am 2 to show up in solar neutrinos. If 
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however, L e is a spontaneously broken quantum number, the associated Nambu-Goldstone 
boson, £ e , will acquire a classic field configuration which may be strong enough to produce 
a back reaction with the net effect of improving the adiabaticity of the v e «-> z/ M transitions. 
The final result is a v e energy spectrum harder than expected. 

In 1987, the existing detectors were only able to detect electron anti-neutrinos but 
the now operating Super-Kamiokande and SNO experiments will be capable of detecting 
supernova v e events. The analysis of the energy distribution can in principle reveal or put 
limits on that kind of effect. 

The scenarios of neutrino mixing change considerably if one considers the evidences 
from atmospheric and terrestrial neutrino experiments (for a review see JT[). The atmo- 
spheric neutrino anomaly and the zenith angle dependence observed by Super-Kamiokande 
can be explained by — > v T oscillations, best fit ||53] Am 2 = 5 x 10~ 3 eV 2 , sin 2 29 = 1. 



The alternative — > v e is excluded by the CHOOZ limits fl26fl . This can still accommo- 
date v e — > or v e — > v T as solar neutrino solutions but that is no longer true if one takes 
in consideration the evidence from the Liquid Scintillation Neutrino Detector (LSND) 
for —>■ v e [27p and — > v e oscillations. The very different Am 2 scales involved 
in LSND (Am^ > 0.2 eV 2 ), atmospheric and solar neutrinos call for a fourth flavor - a 
sterile neutrino. In that picture solar u e s oscillate into the sterile u s . 

We now examine the consequences for supernovae neutrinos always assuming the non- 
adiabatic solar neutrino solution, ignoring for definiteness the possible mixing between v s 
and Up or v T . The oscillation pattern is the following: 1) MSW conversion of v e — > z/ M 
with LSND Am 2 ; 2) sequential oscillation — > v e — > u s , the first a LSND transition, the 
second a solar v process. The outcome is a hard spectrum for v e depleted by v e — > u s , but 
only partially because of the non-adiabaticity of this transition. If alternatively, a NG 
field £ e exists (created by u e — > u^, — > u e and u e — > v a ), it improves the adiabaticity of 
—>■ v s causing a v e depletion stronger than predicted by SM interactions. 

If one repeats the analysis with other scenarios of v mixing the effects will be different 
in detail but with one thing in common: the signature of NG fields is a surprise i.e., an 
oscillation pattern not consistent with the v mixing derived from solar, atmospheric and 
terrestrial v experiments. It should be kept in mind that the situation turns more complex 
and rich if there is mixing between the three NG bosons associated with the three lepton 
flavors, a very natural feature if they are all spontaneously broken. This was explored in 

A point that cannot be overstressed is that the Nambu-Goldstone fields are propor- 
tional to the rate of charge violation processes and therefore to the very reaction rates. 
In the case of neutrino oscillation this manifests as a strong dependence of the Majoron 
fields on the magnitudes of the neutrino luminosities. In the numeric simulation I chose 
10 52 ergs/s for v e and 7 x 10 51 ergs/s for u^, values produced and even exceeded in the 
about half a second that lasts between the neutronisation v e burst and the supernova ex- 
plosion |6|, |29| , |30| . The neutrino luminosities decay afterwards in a time scale of 1 second, 



or rather 4 seconds P0| , pT| , as indicated by SN 1987A events |p!q| . The highest luminos- 
ity happens during the first v e burst - above 10 53 ergs/s in the peak [|6| [2"9| , |30f - and the 



Majoron field may be even stronger then. However, the time scale of the rise and fall 
of the u e signal is about 5 ms, too short to authorize a stationary approximation in the 
calculation of £ e . In fact, distances of the order of 10 10 cm in such a period of time are 
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beyond the light cone and a special study is required. 

The effects of the Majoron fields on the neutrino spectra, if any, will be observed 
in a shorter or longer interval of time depending on the actual scale of lepton number 
symmetry breaking. The observation of such a correlation with the flux magnitudes, by 
itself a signature of the NG fields, would thus provide a measurement of the scale of 
spontaneous symmetry breaking. 
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